Investigation of Heat Conductivity in Relativistic Systems using a Partonic Cascade 
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Motivated by the classical picture ol heat flow we construct a stationary temperature gradient in a 
relativistic microscopic transport model. Employing the relativistic Navier-Stokes ansatz we extract 
the heat conductivity k for a massless Boltzmann gas using only binary collisions with isotropic 
cross sections. We compare the numerical results to analytical expressions from different theories 
and discuss the final results. The directly extracted value for the heat conductivity can be referred 
to as a literature reference within the numerical uncertainties. 



I. INTRODUCTION 

Ultrarelativistic heavy ion collisions (HIC) at the Rela- 
tivistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC) create a hot and dense system of strongly 
interacting nuclear matter that only existed in nature a 
few microseconds after the Big Bang [IH3| . At such ener- 
gies, quarks and giuons are deconfined and a new state of 
matter is formed, the Quark-Gluon Plasma (QGP). The 
main experimental discovery made at RHIC and LHC is 
that this novel state of matter behaves as a strongly cou- 
pled plasma, with the smallest viscosity to entropy den- 
sity ratio, r]/s, ever observed [IHZ]- Currently, relativistic 
dissipative fluid dynamics is the main theory employed 
to describe the space-time evolution of the QGP formed 
in HIC. 

The inclusion of dissipative effects in the description of 
the QGP started only a few years ago. To date the ma- 
jority of studies have focused on investigating the effects 
of the shear viscosity in the time evolution of the QGP 
and in extracting its magnitude from HIC measurements 
[5i-.21, . Nevertheless, there are other sources of dissipa- 
tion that might play a role in the fluid-dynamical de- 
scription of HIC, such as bulk viscous pressure and heat 
flow. While bulk viscous pressure effects have already 
been subject to some amount of investigation p2ti32j . 
the effects of heat flow have been largely ignored. 

The main reason for neglecting heat flow is that most 
fluid-dynamical calculations in the field attempt to de- 
scribe the QGP only at midrapidity and very high ener- 
gies, where baryon number and its corresponding chem- 
ical potential are approximately zero. However, when 
trying to describe the QGP at forward rapidities or at 
smaller collision energies, such as the ones probed in the 
RHIC low energy scan, baryon number can no longer be 
ignored and heat conduction or baryon number diffusion 
might play a more decisive role. On the other hand, be- 
fore implementing heat flow in heavy ion collision simu- 
lations, it is useful to study it with more detail in simpler 
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systems and check how well we are able to describe it in 
such cases. Such types of studies were already started in 
Ref.[33]. 

In order to describe heat flow, one must know at least 
the heat conductivity coefficient. Currently, this is not 
known with the desired precision even for relativistic di- 
lute gases. As a matter of fact, there are several expres- 
sions in the literature for this transport coefficient, e.g. 
from Israel-Stewart theory [34] , resummed transient rel- 
ativistic fluid dynamics |35j , or Chapman-Enskog theory 
|36) . All of these methods give slightly different results 
and it is useful to know which one agrees best with the 
underlying microscopic theory in the dilute gas limit. Re- 
cently, this task has been accomplished for the case of the 
shear viscosity coefficient [571 [55] . 

In this work, we investigate the heat flow of a station- 
ary relativistic dilute gas. Our main purpose is to obtain 
a precise expression for the heat conduction coefficient 
in the kinetic regime. In Ref. |38| . this was done by im- 
posing a stationary velocity gradient as a boundary con- 
dition, and then waiting long enough for the system to 
achieve the Navier-Stokes limit. Once the Navier-Stokes 
limit was obtained, the shear viscosity coefficient was ac- 
curately extracted as the proportionality coefficient re- 
lating the shear-stress tensor to the shear tensor. Here, 
we apply the method used in Ref. [38" , referred to as sta- 
tionary gradients method, to extract the heat conductiv- 
ity of a dilute gas described by the relativistic Boltzmann 
equation. We solve the relativistic Boltzmann equation 
numerically using the transport model BAMPS (Boltz- 
mann Approach for Multi-parton Scatterings) described 
in Refs. [39]. For the sake of simplicity, we shall restrict 
our calculations to a classical gas of massless particles, 
considering only elastic binary collisions with a constant 
isotropic cross-section. 

The paper is organized as follows: In Sec. |Tl]we intro- 
duce the basic definitions of relativistic hydrodynamics. 
In the following Sec. |III| we give an overview of stationary 
gradients. The numerical transport model BAMPS we 
use in this work is introduced in Sec. IIVI Next in Sec. |V| 
we derive an analytical expression for the shape of the 
gradients. In Sec. |VI| we show the method to extract the 
heat conductivity using BAMPS. The results obtained 
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from the numerical calculations are shown in Sec. |VII[ 
where we also discuss the comparison to the analytical 
values. Finally, we give our conclusions in Sec. |VIII| 

Our units are h = c = k = I: the space-time metric is 
given by g''" = diag(l, -1, -1,-1). 



II. BASIC DEFINITIONS 

In relativistic kinetic theory a dilute gas of particles is 
characterized by the invariant particle distribution func- 
tion f{x,p). The macroscopic quantities are obtained 
from the moments of this distribution function. The first 
two moments of f{x,p) correspond to currents of con- 
served quantities: the first moment of f{x,p) leads to 
the particle four-flow (here we consider only binary col- 
lisions and, therefore, particle number is conserved). 
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gd^P aft N 



(1) 



while the second moment corresponds to the energy- 
momentum tensor. 



J^P'P fiP^-)- 



(2) 



Above g is the degeneracy factor, and = {E,p) is the 
four-momentum of the particle. Since we consider the 
masslcss limit, E — \p\. 

The particle four-flow and energy-momentum tensor 
can be decomposed with respect to an arbitrary normal- 
ized time-like four-vector, u'^ = 7(1, v) , where u'^Up, = 1. 
The most general decomposition [36] reads 

N^' ^nut" + V^', (3) 
T'"' = eu^'u•' - PAf"" + W^u" ^Wu^ ^ tt^"^ , (4) 

where A'^'^ — g^^ — u^vy is a projection operator onto 
the 3-space orthogonal to -u^, is the particle diffusion 
four-current, is the energy diffusion four-current, and 
TT^'^ is the shear-stress tensor. Since we are consider- 
ing massless particles, the bulk viscous pressure is zero 
and was omitted from the decomposition above. Further- 
more, we introduced the local rest frame (LRF) particle 
number density 



the LRF energy density 

e EE w^T^'^u^ , 
and the isotropic pressure 

p = -A^,r^73. 



(5) 
(6) 
(7) 



Notice that the particle diffusion four-current can be 
written in terms of as 



(8) 



while the energy diffusion 4-current can be expressed in 
terms of T^^^ as. 







(9) 



The heat flow is defined as the difference between 
energy diffusion and enthalpy diffusion. For the system 
we are considering, it reads 



(10) 



where /i = (e + P)/n is the enthalpy per particle. 

Without loss of generality, we use Eckart's definition 
of the four- velocity [ID], 
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By construction, the Eckart definition of the velocity field 
makes the particle diffusion current vanish, — 0, and, 
consequently, the heat flow in the Eckart frame is exactly 
given by the energy diffusion, i.e. 



(12) 



We remark that the definition of the velocity field will 
not affect the value of the heat conduction coefficient. 

For a gas of classical particles, the thermodynamic 
pressure and particle number density are connected to 
the temperature, T, as follows. 



p = nT. 



(13) 



Since we consider a massless Boltzmann gas, the thermo- 
dynamic pressure p is equal to the the isotropic pressure 
P. The fugacity, A = exp (/i/T), is given by 



A 



(14) 



''cq 



where rtcq — gT^/n'^ is the equilibrium particle number 
density with vanishing chemical potential. 



III. STATIONARY TEMPERATURE 
GRADIENTS 

In Navier-Stokes theory [36^ heat flow is proportional 
to the gradient of the thermal potential, a — fi/T, 



P 



T 



P 



(15) 



where k is defined as the heat conduction coefficient, 
is the ordinary four-derivative, D = u^d^ is the comov- 
ing derivative, and = l^^^dv — — u^D is the 
space-like four-gradient. For the sake of simplicity, we 
assume the system to be homogeneous in the y and z 
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plane (here referred to as the transverse direction) and 
resolve only the dynamics in the cc-direction. We also 
wait long enough so that the system achieves a station- 
ary solution, i.e, dt = 0- In such stationary limit and 
with the spatial symmetries described above, the four- 
derivative simplifies to 9^ = {0,dx,0,0). 

In the stationary limit, conservation of momentum dic- 
tates that the pressure gradient must vanish. Further- 
more, for small velocities, ^ 1, the operator re- 
duces to 



(16) 




order to maintain no pressure gradient, we set the fugac- 
ities to A^^^'' = 1.0, and AJf" = 7.72. The chosen values in 
the reservoirs imply that, p'^f^ — p^^^ — 13.25 GeV/fm'^. 
The collective velocity in both reservoirs is set to zero. 




Then the z-component of the heat flow can be cast in the 
following simple form 



(17) 



Figure 1. A sketch of the setup in BAMPS. In the left and 
right sides (x-direction) of the static box we introduce ther- 
mal reservoirs with fixed chemical potential and temperature. 



In this simplified scenario, the heat conduction coeffi- 
cient can be extracted as the proportionality coefficient 
between the heat flow and the gradient of temperature. 



IV. THE PARTONIC CASCADE BAMPS 

In this work, the relativistic Boltzmann equation is 
solved numerically using the BAMPS simulation, devel- 
oped and previously employed in Rcfs. [38, 39, 41 44]. 
This partonic cascade solves the Boltzmann equation. 



p''d^f{x,p) = C[f] 



(18) 



for on-shell particles using the stochastic interpretation 
of transition rates. In this study we consider only binary 
collisions with constant isotropic cross sections. In order 
to reduce statistical fluctuations in simulations and to 
ensure an accurate solution of the Boltzmann equation 
( 18 1 a testparticle method [39 is introduced: The parti- 



cle number is artificially increased by multiplying it by 
the number of testparticles per real particle, Nt^st- The 
physical results are not affected by this step. 

The simulation of the relativistic Boltzmann equation 
is performed in a static box. The transverse plane (y- 
z plane) of the system is assumed to be homogeneous. 
This is maintained by realizing the collisions of particles 
against the boundaries of the static box as elastic colli- 
sions off a wall. In the x-direction the boundaries of the 
box are fixed to have a constant gradient in temperature. 
In practice, this is implemented in BAMPS by removing 
the particles colliding with the wall in the x-direction. 
Independent of the absorption, the reservoirs emit parti- 
cles with fixed temperature and chemical potential. The 
reservoirs in the left and right boundaries of the box in 
the x-direction, as sketched in Fig[l] are defined to have 
fixed temperatures, T^^ and T^^, and fixed fugacities. 



A^** and A^^'', respectively. Their values in our calcula- 



The set of boundary conditions will lead to a stationary 
temperature gradient in the x-direction. Then the heat 



conduction coefficient can be extracted using Eq. (17 1 



The required time for the system to reach its stationary 
state depends on value of the cross section, the temper- 
ature and chemical potential of the thermal reservoirs, 
and the initial state of the system. 

Similar to the work in [38] , the time the profile needs to 
reach the stationary profile is proportional to the inverse 
mean free path. This implies that the more diffuse the 
system is, the faster it reaches the final stationary profile 
and vice versa. In order to save computational time, we 
always initialize the system as close as possible to its 
stationary state. 

The first step in computing the heat flow is extracting 
the components of the particle four-flow iV^ and energy- 
momentum tensor T'"'. Space in BAMPS is discretized 
in small volume elements, Vc- The distribution function 
f{x,p) of this volume element is reconstructed from the 
momenta distribution of the particles inside it. In this 
scheme, the and T'^'' are computed via the discrete 
summation over all particles within the specific volume 
element and divided by the testparticle number: 



N''{t,x) = 



1 



T^"'{t,x) 



(19) 
(20) 



where Nc is the total number of particles inside the corre- 
sponding volume, t is the time, and x is the space coordi- 
nate (defined to be in the center of the volume element). 
Since we shall be considering only stationary solutions, 
the t dependence of the current will be ommited in the 
following sections. We also mention here, that we aver- 
age ( 19 ) and ( 20 ) over many events in order to reduce 



tions are chosen to T]^^ 



0.5 GeV, Tf'^ 



0.3 GeV. In statistical fluctuations. 
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Using Eqs. (19) and (20) we can determine all neces- 



sary macroscopic quantities. Then, using Eqs. ( 19 ) and 



(11), we compute the four- velocity in the Eckart frame 
and extract the particle number density, n, and the pro- 
jection operator A''". Finally, we compute the heat flow 
current using Eqs. ([9| and (12). 



V. ANALYTICAL EXPRESSION FOR THE 
GRADIENTS 

In this work we aim to extract the gradient of a macro- 
scopic quantity, such as the temperature T or the parti- 
cle density n. The gradients are usually computed using 
finite-difference methods which demands a large amount 
of statistics in order to obtain a sufficiently smooth pro- 
file. Therefore, if one is able to obtain the general form 
of the solution that such gradients must satisfy in the 
stationary regime, it would save a huge amount of com- 
putational runtime. 

As already shown in [38] , quantities that are conserved 
in collisions show a linear behavior between the thermal 
reservoirs. For the boundary conditions implemented in 
this work, this is realized for the particle density. 

For the type of boundary condition considered in this 
work, the particle number density satisfies the following 
solution in the stationary regime 



where 



n{x) 



and 
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-.res 



^rcs 
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L-f 2A 



mfp 



(21) 



are the particle number density in 



the left and right reservoirs, respectively, and L is the 
size of the static box in the x direction. We remark that 
this solution might be modified when inelastic collisions 
are included. 

In general, the above solution does not hold near the 
boundaries (in the x-direction) of the static box, where 
the particle number density is discontinuous (see Ref. 
[38j for details). If the mean free-path, Amfp, of the par- 
ticles is large compared to the size of the box, L, such 
discontinuity will lead to deviations from the above so- 
lution. Here we tackle this problem by making sure that 
L ^ Amfp. This will reduce such finite-size effects to 
the minimum amount possible and guarantees that the 
solution in Eq. (21) provides a very good description of 



the particle number density in most parts of the system. 
Using Eq. ( 21 ) and the fact that the thermodynamic pres- 



sure is constant inside the box, it is straightforward to 
obtain the temperature profile. 



T{x)^p 



-.res 
^R 



^res 



L + 2X 



x + 



-.res 
^R 



-.res 



mfp 



In Fig. [2j we show the particle number density (a), 
thermodynamic pressure (b), and temperature (c) pro- 
files computed with BAMPS (dots) at i = 50 fm and 
(722 = 10 mb, with the boundary conditions specified 
in the previous section. The blue curves correspond to 



> 
O 



45.0 



40.0 



35.0 



30.0 



13.6 



13.4 



13.2 



13.0 



0.50 



0.45 



> 

O 0.40 



0.35 




0.30 



BAMPS 
analytic 



reservoir value 



reservoir value 



-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 
x[fm] 

Figure 2. The numerical solutions extracted from BAMPS for 
particle density n, pressure P and temperature T. The cross 
section is set to (J22 = 10 mb. The solutions are shown at very 
late times, t = 50fm/c, where the profiles are almost static. 
The reservoir values are chosen on the left (right) to T]^'"' = 
0.5 GeV and XT = 1-0 (TF' = 0.3 GeV and \T = 7.72). 
We compare each profile with the analytical results derived 



(22) from Eq. (211. The results are averaged over 500 events. 



the analytical solution described above. It is clear that 
Eq. (21) does in fact reproduce the solutions obtained 



numerically with BAMPS and that finite-size effects are 
negligible. Also, we confirm that the thermodynamic 
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pressure is constant, as it should be. 



VI. EXTRACTION OF THE HEAT 
CONDUCTIVITY 

Since the gradients of the particle number density and 
temperature can be computed analytically from Eqs. ( 21 ) 



and (22), and the pressure p is by construction constant, 



the extraction of the heat conductivity from BAMPS sim- 
ulations becomes straightforward. Using n{x) = ax + b 
for the particle density and the relation ( 13 1 we can sim- 
plify Eq. pTl) to 



q = up 



{ax + by 



(23) 



Finally, the heat conductivity can be calculated using 



(ax + bf 
ap 



(24) 



where the heat flow is directly extracted from BAMPS 
and the constants a and b as well as the constant pressure 
p are analytically known. 

Figure [3] shows the heat conductivity as function of 
a;, at f = 5 fm. We use for the cross section (T22 — 
43 mb. It is clear that k is constant over space, which 
confirms that we reached the asymptotic solution. In 
Fig. |4] we show the heat conductivity as function of time 
t. The system relaxes locally on a very short time scale 
to its stationary solution, and the heat conductivity is 
afterwards basically constant over time. We remark, that 
the relaxation time is proportional to the mean free path, 
Amfp (see Ref. 

In order to extract a precise value for heat conductivity 
K, we take an average over the values obtained in the 
whole system and at all time steps that had reached the 
asymptotic state. 



VII. RESULTS 
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Figure 3. Heat conductivity coefficient extracted from 
BAMPS shown at every point in space. The mean value 
with error bar is also displayed. The results are shown at 
t = 5 fm/c for (722 = 43 mb. 500 events were averaged. 
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Figure 4. Heat conductivity coefficient extracted from 
BAMPS shown at different times averaged over the all po- 
sitions in space. The results are are shown for (T22 = 43 mb. 
500 events were averaged. 



In kinetic theory, there are several different methods to 
compute the transport coefficients appearing in the fluid- 
dynamical equations of motion, e.g. Chapman-Enskog 
theory and the many variations of the method of mo- 
ments. The problem is that each method predicts dif- 
ferent expressions for the transport coefficients. Even 
though these differences are not so large for the shear 
viscosity coefficient, which has been discussed in details 
in Refs. [37J|38j, they can be significantly large for heat 
conductivity, as will be discussed in the following. The 
results from BAMPS can help to clarify which method 
reflects more reliably the underlying microscopic theory. 

The method of moments, initially developed by H. 
Grad for nonrelativistic systems [45] and further ex- 
tended to describe relativistic systems by several au- 
thors, is one of the most commonly employed methods 



in heavy ion collisions. Traditionally, this method is em- 
ployed in the relativistic regime together with the so- 
called 14-moment approximation, originally proposed by 
Israel and Stewart [34J- In this scheme, the momentum 
distribution function, /(x,p), is expanded in momentum 
space around its local equilibrium value in terms of a se- 
ries of Lorentz tensors formed of particle four-momentum 
k'^ , i.e., 1, fc'^, k'^k'^ , .... The 14-moment approximation 
consists in truncating the expansion at second order in 
momentum, i.e., keeping only the tensors 1, fc^, and fc^fc"^ 
in the expansion, leaving 14 unknown expansion coeffi- 
cients. The coefficients of the truncated expansion are 
then uniquely matched to the 14 components of TV and 

Israel and Stewart obtained the equations of motion 



6 



of fluid dynamics and, consequently, the microscopic ex- 
pressions of the transport coefficients appearing in such 
equations, by substituting the momentum distribution 
function truncated according to the 14-moment approxi- 
mation into the second moment of the Boltzmann equa- 
tion, 



/ 



(Pk. 
1" 



By projecting this equation with A"ua, Israel and Stew- 
art obtained an equation of motion for the heat flow and 
a microscopic formula for k [331 SS] • For a massless and 
classical gas of hard spheres, one obtains the following 
expression for the heat conductivity [36l |46] , 



2 

^22 



(25) 



Note that Israel-Stewart's 14-moment approximation 
leads to ambiguous results since it can be substituted in 
any moment of the Boltzmann equation [46] . By chang- 
ing the moment in which the 14-moment approximation 
is replaced, the equations of motion remain with the same 
general form, but the transport coefficients become dif- 
ferent. For example, using the choice of moment applied 
in Ref. |47j . one obtains a quantitatively different result 
for K, 



3 

CT22 



(26) 



Naturally, other choices of moment will lead to even more 
different results, but we shall not list them all here. 

Recently, the derivation of fluid dynamics from the 
method of moments was extended in order to remove the 
ambiguities of the 14-moment approximation |35j . The 
main difference between the 14-moment approximation 
and the theory derived in Ref. |35| , namely Resumnied 
Transient Relativistic Fluid Dynamics (RTRFD), is that 
the latter does not truncate the moment expansion of 
the momentum distribution function. Instead, dynamical 
equations for all its moments are considered and solved 
by separating the slowest microscopic time scale from 
the faster ones. The resulting fluid-dynamical equations 
are truncated according to a systematic power-counting 
scheme using the inverse Reynolds and the Knudsen num- 
bers. The values of the transport coefficients of fluid 
dynamics are obtained by re-summing the contributions 
from all moments of the momentum distribution func- 
tion, similar to what happens in Chapman- Enskog theory 
[48] . In practice, the transport coefficients are given by 
a summation of terms, each term originating from a spe- 
cific moment of f{x,p). Once the transport coefficients 
have converged, the summation can be truncated and the 
moments that are irrelevant can be dropped out. For the 
case of a massless and classical gas of hard spheres, the 
first three moments are enough to obtain a convergent 
value for k, and one obtains 1351 



2.5536 

0'22 



(27) 



Including the next moment would only lead to ~ 1% 
corrections to this expression. 

Finally, heat conductivity can also be computed using 
Chapman-Enskog theory, as was done in Ref. [|36 for a 
gas of classical and massless hard spheres, 



2.44 

^■22 



(28) 



Note that the convergence of the above value was not 
investigated in Ref. ^36j and, consequently, it is not pos- 
sible to infer how precise this result is. 

We found excellent agreement to the result of Ref. [3S] , 
i.e., Eq.(27). This can be seen in Fig. [sj where re- 
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Figure 5. The dimensionless quantity na derived from differ- 
ent fluid dynamical theories (lines) compared to the results 
extracted directly from BAMPS (dots). 



suits obtained using BAMPS are compared to Eqs. (25 1, 
(261, (27) and (28 1. The original Israel-Stewart theory 
(25) yields values roughly ^ 22% too low, while the 14- 



moment approximation with the choice of moment de- 
scribed in Ref. [17| leads to a result that is ^ 17% too 
high. The transport coefficient computed via Chapman- 
Enskog theory is close to the value obtained from micro- 
scopic theory, being ~ 4% below it. However, it is not 
clear whether this can be considered as a failure of the 
Chapman-Enskog theory since the precision of that re- 
sult was not clearly stated and it may be possible that 
it can be improved to match what was extracted from 
the Boltzmann equation. Finally, the heat conductiv- 
ity coefficient computed in Ref. [35 was basically the 
same as the one computed by BAMPS, within the errors 
of both the theoretical calculation (originating from the 
truncation of the moment expansion in the expression 
for k) and the numerical one (originating from the finite 
statistics of the calculation). Nevertheless, all theoreti- 
cal calculations predicted a 1/ 0-22 dependence of the heat 
conductivity coefficient on the cross section; a fact also 
confirmed by BAMPS, see Tab. [ij and Fig.[5| 
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0"22 


K [1/fm^] 


0.043 mb 


609.695698 ± 8.4513 


0.43 mb 


59.588731 ± 0.4694 


4.3 mb 


5.943535 ± 0.2041 


43 m& 


0.607672 ± 0.0411 


430 mb 


0.060449 ± 0.0013 



Table I. Numerical results for the heat conductivity coeffi- 
cient K for various elastic cross sections over several orders of 
magnitude. 

VIII. CONCLUSION AND OUTLOOK 

In this work we extracted the heat conductivity coeffi- 
cient for a dilute gas of massless and classical particles de- 
scribed by the relativistic Boltzmann equation. For this 
purpose we employed the microscopic transport model 
BAMPS in a static system using only binary collisions 
and a constant isotropic cross section. For this setup, we 
established a stationary temperature gradient using ther- 
mal reservoirs. In order to simplify the calculations we 
derived an analytical expression of the expected profile 
to obtain the gradients of the temperature, meanwhile 
the heat flow was extracted directly from BAMPS by 
the general decomposition of the particle four-flow and 
energy momentum tensor. Using the relativistic Navier- 
Stokes theory, which is valid for a stationary system and 
small gradients, we extracted the heat conductivity to a 
very high precision. We then compared this result with 
several theoretical predictions for this transport coeffi- 
cient, each originating from a different derivation of fluid 
dynamics from the underlying microscopic theory. 

The numerical simulations of the Boltzmann equation, 
realized in this paper with BAMPS, were able to distin- 
guish between all these different theoretical results and 
clearly point out which one was in better agreement with 
the underlying microscopic theory. While Israel- Stewart 



theory, i.e., the 14-moment approximation, performed 
rather poorly in the description of heat flow, new ex- 
tensions of the method of moments, i.e. RTRFD, were 
able to provide an improved description of this transport 
coefBcient, showing a very good agreement with the co- 
efficient extracted from BAMPS. The heat conductivity 
computed using Chapman-Enskog theory was not able 
to precisely describe the simulation from BAMPS, al- 
though it was not as far off as the predictions using the 
14-moment approximation. Nevertheless, this disagree- 
ment might be originating from a poor implementation 
of the Chapman-Enskog theory, which was not properly 
checked for the convergence of the transport coefficient. 
The extracted value for the heat conductivity 

2.59 ±0.07 

'^BAMPS = • (29) 

can be referred to as a literature value, within the numer- 
ical uncertainties, for the simple case of binary collisions 
with an isotropic cross section. It remains as a future 
task to extend this work to include inelastic scatterings, 
where particle production and annihilation have to be 
taken into account. 
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